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^--^ Aims. We present photoionization calculations for the spatially extended absorbers observed in front of the extended emission 
\^ line spectrum of two high redshift radio-galaxies: 0943-242 {z e — 2.922) and 0200+015 (z e = 2.230), with the aim of reproducing 
f****) , the absorber column ratio, Now /Nni- 

Methods. We explore the effects of using different UV continua in the photoionization calculations. A comparison is drawn 
between the absorber in 0200+015 and the two absorbers observed near the lensed Lynx Arc Nebula at redshift 3.36, which 
present very similar Nciv/Nm ratios. 
Q ' Results. We find that hot stars from a powerful starburst, or a metagalactic background radiation (mbr) in which stars dominate 
S—j \ over quasars, are equally successful in reproducing the observed Ncrv /Nm, assuming subsolar gas metallicities for each absorber. 
■ These softer SEDs eliminate the difference of a factor 1000 in metallicity between the two absorbers encountered in earlier work 
where a powerlaw SED was assumed. The detection of continuum flux in 0943-242 suggests that the level of ionizing photons is 
consistent with a stellar ionizing source. 

Conclusions. If the mbr is responsible for the ionization of the radiogalaxy absorbing shells, their radii (if spherical) would be 
large (> 100 kpc) and their mass would be huge > 10 12 Mq, implying that the feedback mechanism initiated by the central 
galaxy has caused the expulsion of more baryonic mass than that left in the radiogalaxy. If, as we believe is more likely, stellar 
ionizing sources within the radio galaxy are responsible for the absorber's ionization, smaller radii of ~ 25 kpc and much smaller 
masses (~ 10 s — 10 10 Mq) are inferred. This radius is consistent with the observed transition in radio source size between the 
smaller sources in which strong Hi absorption is almost ubiquitous, and the larger sources where it is mostly lacking. Finally, 
we outline further absorption line diagnostics which could be used to constrain further the properties of the haloes and their 
source of ionization. 
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1. Introduction 

A prominent characteristic of high-redshift radio-galaxies 
(hzrgs) at z > 2 is their spatially extended line emission 
regions (hereafter eelr), which are often luminous in Lya 
(> 10 44 ergs -1 ) and extended over several to tens of kpc. 
The excitation mechanism for the emission gas is either 
shock excitation by jet material or AGN photoionization 
(the presence of NvA1240 line emission precludes stellar 
photoionization). The eelr is kincmatically active, with 
FWHM reaching 1000 km s~ 4 . With observations of a sam- 
ple of ffiRGs, VanOjik etal. (1997;V097) discovered that, 



when observed at intermediate resolution (1-2 A), the ma- 
jority of HzRGs with small radio source sizes (< 50 kpc) 
exhibit narrow Lya H I absorption. This absorption is su- 
perimposed upon the Lya emission with a spatial extent 
comparable to that of the eelr. In addition to Lya, the 
CivAA1549 doublet has also been observed in absorption 
in two hzrgs, superimposed on the Civ emission line, 
firstly in 0943-242 (z e = 2.922) (Binette etal. 2000, here- 
after BOO) and secondly in 0200+015 (z e = 2.230) (Jarvis 
etal. 2003, hereafter J03). Building on the results of BOO 
and J03, in the present paper we examine the excitation 
mechanism of the large scale absorbing haloes in greater 
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detail by exploring photoionization by a variety of differ- 
ent spectral energy distributions (hereafter sed). 

The basic structure of the paper is as follows. In the re- 
mainder of section 1 we review our current understanding 
of hzrg absorbers, focussing on the distribution, ioniza- 
tion and metallicity of the absorbing gas and the specific 
problems which motivate our current study; an insight- 
ful comparison is made with the absorbers in the Lynx 
Arc Nebula (lan), a gravitationally-lensed Hn galaxy at 
z = 3.357. In section 2 we summarise the observational re- 
sults we aim to reproduce, namely the Ncrv/Nm ratio in 
the aforementioned hzrgs and the LAN. Section 3 describes 
the mappings ic code and our assumptions concerning the 
photoionizing SEDs. Section 4 presents the results of these 
calculations and in section 5 we assess their implications 
for the origins of the absorbers and their compatibility 
with other observables. Finally, in section 6 we present 
some additional absorption line diagnostics which may in 
future help to discrimate between the proposed scenarios. 

1.1. Shell-like structure for the hzrg haloes 

Among the hzrgs with small radio sources (< 50kpc), the 
detection rate of associated absorption systems is 90% (9 
out of 10 hzrgs in the V097 study) while it is only 25% 
for larger radio-sizes. The fact that the absorption extends 
over the whole background eelr emission favours a shell- 
like geometry for the absorption systems rather than a 
conglomerate of individual clouds, as proposed initially by 
V097. In Sect. ll.3l we give further indications as to why we 
retain the simplifying assumption of a simple shell struc- 
ture in the current work. Because the density per unit red- 
shift of the strong absorbers (iVgi > 10 18 cm -2 ) around 
HzRGs was found to be much higher than that given by 
the statistics of intergalactic medium (IGM) absorbers at 
large, V097 inferred that they belong to the environment 
of the parent HzRG rather than to the IGM. The density 
of the thinnest absorbers (< 10 15 cm -2 ) around hzrgs on 
the other hand is comparable to that of Lya forest ab- 
sorbers in the IGM, as more recently shown by Wilman 
etal. (2004: W04). It is conceivable that the physical con- 
ditions in the thin 1 HzRG absorbers are indistinguishable 
from those operating within typical IGM Lya forest ab- 
sorbers. The available data, however, are still insufficient 
to confirm or refute this proposition. 

The rarity of absorbers among hzrgs with radio-sizes 
larger 2 than 50kpc suggests that the typical lateral di- 
mensions of the shell (in the plane of the sky) might be 
^ 50kpc. The proposed interpretation is that, as the AGN 
jet expands beyond such size, the bow-shocks overtake the 
shells and disrupt them. This is the first scenario, which 
we label A or "inner shell scenario" . If valid, it suggests 

By "thin" we refer to a small Hi column density, not to a 
small physical size. 

2 V097 (p. 369) showed that this finding is not the result of 
a selection effect due to the fact that larger radio sources have 
a narrower Lya emission profile. 



that the expansion of the AGN jet cocoon is not the mech- 
anism by which the shells are formed, but rather by which 
they are destroyed. Scenario A favours a shell formation 
mechanism that relies on large-scale outflows generated by 
episodes of massive star formation. Using high dispersion 
data from VLT-UVES, W04 proposed that the absorbers 
in HzRGs probably lie within the core of young galactic 
proto clusters, consistent with observations of their envi- 
ronments (e.g. Venemans et al. 2005; Overzier et al. 2006), 
and may be a byproduct of massive galaxy formation. 
Krause (2005) published hydrodynamical simulations of 
the formation of a shell due to the expansion of a stellar 
wind bow-shock. At a later stage in his model, an AGN 
jet is launched and a jet cocoon builds up. Once the jet 
has extended beyond the initial bow-shock, the jet cocoon 
destroys the shell as it overtakes it. An estimate of the 
timescale for this to occur can be obtained if one follows 
the reasoning of J03, where the radio-size represents a 
kind of internal clock (see Sect. ll.5jl . which characterizes 
not only the radio jet's age but also that of the starburst 
superwind that generates the shells. 

A second possibility is that the rarity of shells among 
HzRGs with large radio-sizes may reflect an older phase in 
which the shells have expanded further out and thinned 
out considerably. This process would eventually render 
them undetectable (using the V097 detection technique) 
when their A/hi columns drops below ^ 10 13 cm~ 2 . This 
is the second scenario, which we label B or "aging shell 
scenario" . In this case, the distance between the shell and 
the parent hzrg is unknown and can be much larger than 
the upper limit size implied by scenario A, as will be dis- 
cussed in 15.11 Scenario B leaves the possibility open that 
some of the shells may result from the expansion of a jet 
bow-shock (e.g. Krause 2002), although the most likely 
formation mechanism of the shells remains stellar super- 
wind, as in scenario A. 

The large scale HzRG absorbers might be analogous to 
the absorbers detected within 3 20-50/i _1 kpc of high red- 
shift galaxies, by Adelberger etal. (2005) using nearby- 
field spectroscopy of background QSOs or galaxies. The 
advantage of hzrg absorber studies is that the intrinsic 
shell outflow velocity is more readily available from obser- 
vations, but not their distance from the parent HzRG (the 
reverse applies to the technique used by Adelberger et al. , 
assuming spherical expansion of the absorbers). 

The HzRG shells share many similarities with at least 
some of Lya-emitting 'blobs' (hereafter LAB; Steidel et al. 
2000) associated with Lyman Break Galaxies (Pettini 
etal. 2001), and which are characterized by an eelr that 
can reach large sizes of up to ~ 100 kpc. An important 
difference is that the radio luminosities are much fainter 
or undetected in the latter case. Using integral field spec- 
troscopy, Wilman et al. (2005) observed the Lya-emitting 
'blob' LAB-2 in the SSA22 protocluster at z = 3.09; 
they discovered a foreground absorber (Ahi — 10 19 cm~ 2 ) 
with remarkable velocity coherence over a projected size of 

3 Where h = ^/lOOkms^ 1 Mpc" 1 . 
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~ 76 x 26kpc. Their interpretation is that a galaxy-wide 
superwind has swept up ~ 10 11 Mg of diffuse material 
from the IGM over a few 10 8 yr. This is a manifestation 
of the 'feedback' mechanism thought to be regulating the 
formation of galaxies. 

1.2. What is ionizing the large scale haloes? 

At least a fraction of the known haloes appear to be highly 
ionized. In both hzrgs in which a C iv AA1549 doublet has 
been observed in absorption (0943-242 and 0200+015), 
the absorption rcdshift corresponds to one of the Lya ab- 
sorbers. These two absorbing haloes therefore contain ion- 
ization species up to at least C +3 . BOO and J03 assumed 
that the Hi and Civ absorption species occur within a 
physically contiguous structure, an aspect discussed fur- 
ther in Sect. II. 31 

The possibility that the hzrg absorbing haloes are pho- 
toionizcd by the hidden nuclear radiation can be ruled 
out. Firstly, because there no observed continuum of suffi- 
cient strength underlying the eelr. This is as expected in 
the quasar-radio-galaxy unification picture (Barthel 1989; 
Antonucci 1993; Haas etal. 2005), in which the nuclear 
ionizing radiation is collimated along two ionization cones, 
which in radio-galaxies lie along the plane of the sky, 
and is therefore invisible to the observer and presum- 
ably to the intervening absorbers as well, unless rather 
contrived gas geometries are postulated. Secondly, BOO 
showed that the Civ/Lya emission- and absorption-line 
ratios in 0943-242 could not be reconciled with any model 
in which the absorption- and emission gas are co-spatial. 
They concluded that the absorbing gas is of much lower 
metallicity and located further away from the host galaxy 
than the eelr. Although they favoured the idea that the 
diffuse metagalactic background radiation (hereafter MBR) 
rather than the parent AGN was responsible for ionizing 
the absorbing haloes, in their calculations BOO and J03 
used a simple powerlaw as a crude approximation of the 
mbr energy distribution. In this paper, we will assume 
more realistic SEDs that take into account the cumulative 
opacity of IGM Lyman limit systems and Lya forest ab- 
sorbers. 

As for the possibility of collisional ionization of the 
shells, we indicated in J03 that this mechanism was un- 
likely in the case of 0200+015, and that steady-state pho- 
toionizing shocks (Dopita & Sutherland 1996) resulted in 
rather large Ahi columns (~ 10 19 cm -2 ), incompatible 
with the small value characterizing 0200+015. In the case 
of 0943-242, the near solar metallicity models of Dopita & 
Sutherland (1996) do not attain the observed AVjrv value 
for shock velocities below 400 km s _1 and above this veloc- 
ity the Ahi column becomes excessive, requiring the shock 
structure to be truncated. As for the photoionized precur- 
sor nebula upstream from the shocks, the SED generated 
downstream by fast shocks is as hard as a powerlaw of 
index a ~ —0.5 up to <; 500 eV (Binette, Dopita & Tuohy 
1985). Therefore, photoionization calculations with a pow- 



erlaw as presented in Sect. l4.2| capture the main features of 
such a precursor. In essence, any hard SED requires super- 
solar metallicities in order to fit the column ratios found in 
0200+015. Finally, calculations to represent the case of a 
collisionally ionized gas slab at temperature T has been ex- 
plored by J03. They find that for 0200+015, the N cw /N m 
column ratio could be reproduced by using roughly solar 
metallicities, provided that T is finetuned to lie around 
10 5 K. Apart from the fact that this metallicity is rather 
large for the redshift considered, it would be difficult to 
explain how the plasma could be maintained at a tempera- 
ture approaching the peak of its cooling curve. This would 
most likely require a yet unknown heating mechanism. 

1.3. The case for a simple scattering screen 

In a morphologically and kinematically complex eelr, 
one cannot readily disentangle photon destruction due to 
line-of-sight absorption from the effects of transmission 
by multiple scatterings. Nevertheless, for the large scale 
absorbing haloes in 0200+015 and 0943-242 (or other 
hzrgs studied by V097), there is no evidence that the 
absorbers share the complexity of the eelr. The results 
suggest that a uniform foreground scattering screen pro- 
vides an adequate description of the "absorbing" haloes in 
hzrgs. Strong evidence for this was provided by observa- 
tions at much higher spectral resolution using VLT-UVES 
(e.g. J03 and W04). In particular, J03 finds that the main 
absorber in 0943-242 remains as a single system of col- 
umn density ~ 10 19 cm -2 over the full size of the eelr, 
being completely black at its base, with no evidence for 
substructure or a multiphase environment. The absorp- 
tion trough is blueshifted by 265 km s -1 with respect to 
the centroid of the background emission profile. This spa- 
tial and kincmatical coherence of the absorber contrasts 
with the chaotic multiphase medium encountered in the 
Galactic ISM or the eelr of hzrgs. It also suggests that 
the absorber is physically separate from the background 
eelr and that it is therefore simply acting as a scatter- 
ing surface, as argued in J03. This clean separation be- 
tween eelr and the absorber simplifies the modelling task 
and justifies the ionization-stratified slab approximation 
adopted in Sect.|H 

1.4. Comparison with the Lynx Arc Nebula: LAN 

Following an independent study of the lensed Lynx Arc 
Nebula 4 at z = 3.357 by some of us (Fosbury et al. 2003), 
it was observed that the column ratios Ncw/Nm in the 
LAN and 0200+015 are very similar. In the calculations 
which follow, we therefore explicitly compare the LAN and 
the HzRGs absorbers, making use of the following insights 
which place the physical conditions in the LAN on a firm 
footing. Firstly, the LAN is an active star-forming object, 

4 The Lynx Arc Nebula is a high-redshift metal-poor gravita- 
tionally lensed H II galaxy that was discovered serendipitously 
by Holden etal. (2001). 
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so we may reasonably assume that the subsolar metal- 
licity that characterizes the emission gas, ~ 10% (fol- 
lowing the work of VM04), also applies to the absorbing 
gas. Secondly, the LAN presents a relatively high excita- 
tion emission line UV spectrum, which photoionization 
by hot stars can reproduce successfully. Photoionization 
by a straight powerlaw 5 , on the other hand, would re- 
sult in the emission of a detectable NvA1240 line (com- 
parable in strength to Niv] A1485 line, see BG03), which 
is not observed. Hence, for this object, we know the ab- 
sorber metallicity and excitation source with some confi- 
dence. Therefore, the successful reproduction of the LAN 
column ratios in Sect . 14.41 using subsolar metallicities and 
photoionization by hot stars, prompts us to consider that 
such an SED might also apply to the hzrg absorbers. 

1.5. Metallicity evolution vs. softer ionizing SED 

With the VLT-UVES, J03 obtained superb spectra of the 
aforementioned hzrgs at ten times the resolution used by 
V097. The spectra confirmed that the main absorber in 
0943-242 exhibits no additional substructure to that re- 
ported by V097, as already discussed. In contrast, a very 
different view of 0200+015 emerges: the single absorber 
with HI column density ~ 10 19 cm -2 seen at low reso- 
lution now splits into two ~ 10 146 cm -2 systems; these 
extend by more than 15 kpc to obscure additional Lya 
emission coincident with a radio lobe. Additional but frag- 
mented absorbers are seen on the red wing of the emission 
line at this position. We recall that gas metallicities as 
high as ~ 10 Zq are required to reproduce the Aciv/Ahi 
ratio in 0200+015 (Sect. 11.51 J03) assuming photoioniza- 
tion by a straight powerlaw. This suggests that the ab- 
sorbing gas has undergone very substantial metal enrich- 
ment. Based on the smaller radio source size in 0943-242 
(26 kpc versus 43 kpc for 0200+015), J03 conjectured that 
the radio source age (as inferred from its linear size) is 
the parameter controlling the evolution of: (i) the struc- 
ture/kinematics of the absorbing halo, through interac- 
tion and shredding of the initially quiescent shells; and 
(ii) its metallicity built-up, through enrichment by star- 
burst superwind triggered concurrently with the nuclear 
radio source. 

Although this age and enrichment scenario (B) re- 
mains an appealing possibility, the large metallicity gap 
inferred by J03 of three orders of magnitude between 
0943-242 (~ 0.01 Zq) and 0200+015 (~ 10 Z Q ) is a cause 
for concern. Here we revisit the issue by exploring alter- 
native ionizing SEDs that would require only a factor of 
ten metallicity enhancement with respect to 0943—242. 
We focus on the case of SEDs from hot stars and the dif- 



5 The possibility that the LAN corresponded to photoioniza- 
tion by a quasar powerlaw that is partly absorbed has been 
investigated by Binette etal. (2003: BG03). In order to work, 
such an absorbed powerlaw model requires a very specific and 
unlikely fine tuning of the parameters defining the eelr and 
the hidden filtering screen. 



fuse mbr with the aim of generating a grid of models for 
comparision with future observations. 

2. The observational dataset 

In this section, we gather together the principal observa- 
tional results which we aim to reproduce, namely the H I 
and C iv column densities for the absorbers in the two 
HzRGs and the LAN. 

The C IV and the Lya absorption columns in 0943-242 
(z c = 2.922) and 0200+015 (z c = 2.230) have been mea- 
sured by various authors (Rottgering et al. 1995; BOO, J03, 
W04). We adopt the values of J03, which are based on 
VLT-UVES observations of both hzrgs. In 0943-242 the 
dominant large scale absorber is characterized by an JVhi 
column of 10 191 cm -2 , which puts it among the group 
of larger Hi columns (see W04). However, the four Lya 
absorbers observed in 0200+015 are rather thin, with 
columns of order 10 147 cm~ 2 . These all belong to the 
group of smaller H I columns haloes, which are much more 
numerous (see W04). In 0943-242, the CivAA1548, 1551 
doublet is observed in absorption at the same redshift 
as the dominant Hi absorber (J03; BOO; Rottgering & 
A4iley 1997) and corresponds to a column of 10 14 6 cm -2 . 
In the case of 0200+015, only one Hi absorber, with 
TVhi = 10 14 7 cm~ 2 shows a corresponding C iv doublet 
in absorption, with AT C rv = 10 14 - 6 cm -2 . The N C w/N m 
column ratios for 0943-242 and 0200+015 are 10~ 4 - 5 and 
10~ 07 , respectively. 

As for the LAN, the two local absorption systems in the 
Lynx Arc have been labelled al and a2 by Fosbury et al. 
(2003) who determined the Hi column to be 1.05 x 10 15 
and 0.60 x 10 15 cm -2 , respectively, and the Civ columns 
to be 0.83 x 10 15 and 1.02 x 10 15 cuT 2 , respectively. 
The iVciv/A^Hi column ratios for al and a2 are therefore 
iQ-o.io and 10 23 , respectively. The similarity of al with 
0200+015 is noteworthy. 

3. Photoionization models and ionizing energy 
distributions 

3.1. mappings ic and gas abundances 

To compute the Nciv/Nm ratio, we have used the code 
MAPPINGS ic (Binette, Dopita & Tuohy 1985; Ferruit et al 
1997). To represent solar abundances, we adopt the set 
of Anders & Grevesse (1989). When varying metallicities, 
we multiply the solar abundances of all elements heavier 
than He by a constant, which we label the gas metallic- 
ity (in units of Zq). For the HzRG absorber, we assume 
a slab geometry illuminated on one side. For each ioniz- 
ing SED that we considered, we calculated the equilibrium 
ionization state of the gas and integrated the ionization 
structure inward until a preselected target value of the 
column Ahi was reached. For the range of parameters 
explored in this paper - where the aim is to reproduce 
the observed JVhi and the Aciv/Ahi column ratios - all 
models of the absorbers turn out to be matter-bounded 
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(0200+015) or marginally optically thick to the ionizing 
radiation (0943-242). 

3.2. Continuum softness 

For a given input SED, the calculations are repeated for 
different values of the ionization parameter 6 in order to 
build a sequence of models in U, starting at the minimum 
value of 0.001. 

It is customary to define the SED's softness using the 
parameter 77, which is the column ratio of singly ionized 
He to neutral H, Nn e n/ Nm- This ratio does not, how- 
ever, uniquely define the SED, as 77 also depends on the 
slab thickness and on U, and not just on the continuum's 
shape (see for instance Appendix A of Fardal etal. 1998). 
In the case of stellar SEDs, 77 varies rather abruptly with 
T e g. For instance, 77 is 1480 for a 71 000 K star while for 
a 80 000 K star its value is only 95. Next to each SED in 
Fig.Q] we indicate the value of r\ calculated between brack- 
ets, assuming Ahi = 10 14 ' 8 cm -2 and U — 0.1. We now 
review the various SEDs displayed in Fig^and used in the 
calculations reported in Sec. El 



3.3. AGN powerlaw SED 

In the case of direct photoionization by an AGN, we as- 
sume a simple power law of index a = —1.0 as in BOO 
(with J v oc v +a ). This SED is represented by the gray 
long-dashed line in Fig.^. 



3.4. sed of the diffuse metagalactic radiation (mbr) 

The integrated ultraviolet flux arising from distributed 
QSOs and/or from hot massive stars (in metal-producing 
young galaxies) is believed to be responsible for maintain- 
ing the intergalactic diffuse gas, the Lya forest and Lyman 
limit systems in a highly ionized state. The spectrum and 
intensity of the diffuse mbr is affected not only by photo- 
electric absorption from intergalactic matter, but also by 
the re-emission from radiative recombinations within the 
absorbing gas itself. In short, QSO absorption- line systems 
are sources, not just sinks of ionizing photons, as shown 
by Haardt & Madau (1996). 

Detailed calculations of the propagation of QSO and 
stellar ionizing radiation through the intergalactic space 
have been presented by Fardal etal. (1998; FGS in the 
figures or footnotes) and the resultant SEDs relevant to 
the current work 7 are shown in Fig.^i. On the one hand, 
we have the metagalactic SED in which only quasars are 



We use the customary definition of the ionization parame- 
ter U = (fin/enn as the ratio between the density of ionizing 
photons impinging on the slab <pn/c and the total H density 
at the face of the slab tih- 

7 The SEDs calculated by FGS are softer than those of Haardt 
& Madau (1996). In their AppendixA, FGS justifies this dif- 
ference on account of a more detailed treatment of the cloud 
opacity and re-emission. 
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Fig. 1. The spectral energy distribution of various ioniz- 
ing sources (see Sect.[3J) as a function of photon energy. 
Panel a: silver long-dashed line: spectral energy distribu- 
tion corresponding to an AGN powerlaw, short dashed- 
linc: diffuse mbr energy distribution from FGS at z a = 2 
comprising only quasars (Q) as sources, continuous line: 
(original) MBR energy distribution from FGS at z a = 3 
comprising both quasar and stellar sources (Q+"*"), dotted 
line: same as solid line except that the flux beyond 54.4 eV 
has been reduced by a factor two. Panel b: continuous line: 
SED of a zero-age metal-free star of effective temperature 
80 000 K, dotted-line: SED of a zero-age metal-free star of 
effective temperature 88 000 K, long-dashed line: SED from 
an evolutionary model of CMK04 corresponding to an age 
of 3.4 Myr and metallicity of 20% solar. Values of 77 are 
given between brackets for each SED fSec. !3.2|l . 



contributing (short dashed line) corresponding to model 
Q2 in their Fig. 7 and, on the other, the SED in which hot 
stars from star forming regions are included (continuous 
line), a model shown in their Fig. 6. In this model, the 
stars are contributing twice the flux of quasars at 13.6 eV. 
The dotted line is a similar SED, except that the flux be- 
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yond 4Ry has been divided by two. It is an ad hoc model 
representing the case in which stars are contributing pro- 
portionally more with respect to quasars (a similar SED 
was also considered by Telfer etal. 2002). 

The sharp drop of flux at 54.4 eV is a characteristic 
of all metagalactic radiation models and is due to the cu- 
mulative opacity of Hell within the IGM. As the IGM 
SEDs extend into the soft X-rays, it is important to in- 
clude the harder radiation beyond ^ 200 eV, otherwise 
the calculated C iv columns are affected, especially when 
U is large. 

If we turn to the values of the softness parameter 
(Sect. 13.21) observed among IGM absorbers, there is a sub- 
stantial dispersion in the values measured by Kriss et al. 
(2001), with 1 < i] <; 1000, which suggests that for a 
fraction of absorbers, stellar ionizing sources might be 
contributing. The possibility of a rather inhomogeneous 
distribution of the SED hardness according to location is 
favoured by the independent study of Smette et al. (2002), 
who find that 20 <> 7} 5000. 

3.5. Stellar ionizing seds 

In the stellar ionizing case (panel b in in Fig.Q, we consid- 
ered metal-free stellar SEDs that approximate those stud- 
ied by Schaerer (2002) with T c g among one of the follow- 
ing values: 42 000, 57000, 71000, 80 000 and 88 000 K. In 
Fig-Ub, we illustrate the cases of the 80 000 and 88 000 K 
seds (the continuous and dotted lines, respectively). As 
in BG03, who presented various photoionization models 
for the lan, we approximate the selected stellar seds, 
using a technique that reproduces the ionizing photon 
luminosities Q(H), <5(He°) and Q(He + ) of the selected 
logT e ff model listed in Table 3 of Schaerer (2002). In a 
similar fashion to Shields & Searle (1978), we derive the 
monochromatic temperatures at the edge boundaries T^ , 
T~ c+ and + ) and then interpolate linearly in log T v 
for all the wavelengths used in the code MAPPINGS ic. We 
equated T^ to T e g and neglected the very small He edge 
present in these atmospheres. This simplified representa- 
tion of a stellar atmosphere provides enough accuracy to 
compute the essential properties of the emission line spec- 
trum. 

We additionally considered an SED derived from the 
stellar evolutionary model of Cervino, Mas-Hesse & Kunth 
(2004, hereafter CMK04), which was used by Villar- 
Martfn etal. (2004; VM04) in their photoionization cal- 
culations of the LAN. The selected SED corresponds to a 
metallicity Z* = 0.20 Zq and an age of 3.4 Myr. It is rep- 
resented by the long-dashed line in Fig.^h We included 
the weak X-ray flux that results from the conversion of the 
kinetic energy of the supernova remnants into X-ray emis- 
sion (it did not have any effect on the results). The stellar 
cluster at that particular age harbours an important popu- 
lation of WR stars and, as shown by VM04, the resulting 
ionizing continuum is sufficiently hard to reproduce the 
emission line strength of the He II A1640 line observed in 



the LAN spectrum. The CMK04 evolutionary models are 
characterized by a powerlaw initial mass function with a 
Salpcter IMF and stellar masses comprised in the range 
2-120 Mq. 

4. Model results 

In this section, we present a grid of photoionization cal- 
culations for comparison with the observed Ncw/Nm ra- 
tios in the two HzRGs and the LAN. In Sect. l4.l| we outline 
our investigative procedure and the format we adopt to 
display the results. Thereafter, we explore the effects of 
using different seds and varying some of the input pa- 
rameters, as follows: (a) powerlaw photoionization is first 
studied in Sect . 14.21 assuming different metallicities; (b) 
in Sect. 14.31 we study various MBR energy distributions in 
which quasars and stars are contributing in different pro- 
portions; (c) in Sect. l4.4l we explore stellar photoionization 
by metal-free atmospheres of varying T e s and by a stellar 
cluster SED containing WR stars. 

4.1. Aims and modelling procedure 

There is a gap of more than four orders of magnitude 
in the N C w/N m ratio between 0943-242 and 0200+015. 
Rather than explain this with a factor ~ 1000 difference in 
absorber metallicity between the two HzRGs as in J03, we 
instead explore alternative SEDs. As stated in Sec. II .51 our 
practical goal is to find an SED which reduces the metallic- 
ity gap to ~ 10 (that is, obtaining a successful model that 
use abundances as low as ~ 10% solar). We do not aim at 
obtaining exact fits of this ratio in each case, but rather 
at establishing an order of magnitude agreement between 
the models and the separate observations of the thin and 
thick absorber categories. For this reason, we only con- 
sider the following four widely spaced metallicities for the 
haloes: 10, 1, 0.1 and 0.01 Zq. The higher the metallicity, 
the larger the Ncw/Nm ratio. The proportionality is lin- 
ear except in the high metallicity regime, where the slab 
temperature structure is somewhat altered, this effect be- 
ing more important in the case of the thick absorbers. No 
attempt is made in this paper to model the background 
eelr spectrum. The metallicity of the eelr gas is much 
higher than (and unrelated to) the absorber's, as discussed 
by BOO. 

The ionization parameter characterizing the models is 
plotted on the abscissa in all figures. It is a free parame- 
ter that cannot be adequately constrained with the limited 
data at hand. The target Nqw/Nhi ratio (y-axis) for each 
observational datum is represented by an horizontal (thick 
broken) line, since U is not known. Our aim will be to find 
models that either cross this observational line or come 
close to it. We consider it unlikely that U is smaller than 
0.001, since C +3 would then be reduced to a trace species. 
It is plausible that it takes on much larger values instead, 
especially in the case of 0200+015 or the LAN, since large 
values of U usually result in larger Aciv/^hi ratios, a 
characteristic of these thinner absorbers. We adopt the 
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1 - 3 0.01 0.1 1 

ionization parameter U 

Fig. 2. The column ratio Nciv/Nm derived from pho- 
toionization by a powerlaw of index —1.0, as a function of 
the ionization parameter U. Each model sequence (thin 
lines) along which U varies represents a slab of fixed 
Nfti column (shown between brackets). Thin black lines 
connect models of constant column Nm = 10 148 cm~ 2 , 
while thin gray lines connect models of constant column 
Nm = 10 19 cm -2 . The gas metallicity is shown using la- 
bels, in units of Zq. The thick horizontal broken lines rep- 
resent four measurements of Aciv/-^hi: the extended ab- 
sorbers in the high-z radio galaxies 0943-242 (in gray) and 
0200+015 (in black) and the two absorbers found in the 
lensed Lynx Arc Nebula [lan] (in black) and labelled al 
and a2 by Fosbury etal. (2003). The ratios for 0200+015 
and al are very similar and have been combined into a sin- 
gle entry. The object's name and the Nm column appear 
to the left and right, respectively, of the corresponding 
horizontal broken line. 

conservative view that most of the difference between the 
thin and thick absorbers may be accounted for by differ- 
ences in the gas excitation, that is, in U rather than by 
metallicity differences only. Future observations of other 
resonance lines might be used to test this (see Sec. 16.1(1 . 

In what follows, the slab models presented are either 
characterized by a small Hi column of 10 148 cm~ 2 as in 
0200+015 and the LAN, or a larger column of 10 19 cm -2 
as in 0943-242. We emphasize that in all the figures (J2HEJl 
of Section^] the black line models only apply to the thin 
absorbers shown at the top while the gray-line models only 
apply to 0943-242 (shown at the bottom). 

4.2. Powerlaw photoionization 

We present photoionization calculations in Fig.[21 for the 
case of an AGN powerlaw of index —1.0. Using a mod- 



erately different index would not significantly alter the 
conclusions reached below. For instance, a steeper index 
— 1.4 would increase the column ratio only by a factor of 
£ 2. 

In the case of the 0943-242 absorber, the models in 
Fig. |21 (thin gray lines) favour abundances much lower than 
solar, that is of order 1% solar. BOO favoured a metallicity 
value of ~ 0.02 Zq. A much lower (higher) ionization pa- 
rameter is a possibility that cannot be ruled out and the 
metallicity would then be higher (lower) than the values 
we considered. 

In the case of the 0200+015 absorber, very high metal- 
licities are favoured by the powerlaw SED, as found by 
J03. This is shown by the thin black lines in Fig.0 which 
suggest a gas metallicity of about ten times solar. We ex- 
pressed concerns about such high values in Sect. 11.51 

We reject the powerlaw SED on account of geometrical 
considerations presented in Sect. ll.2l that led us to rule out 
direct ionization by the nuclear radiation from the AGN. 
The main purpose in reporting powerlaw calculations is 
to provide a convenient comparison with the softer con- 
tinuum shapes explored below. We note how different the 
behaviour of Aciv /Nm is between the thin and thick ab- 
sorber case in Fig. (compare the two long-dashed models 
of solar metallicity). 

4.3. Photoionization by the diffuse mbr 

4.3.1. mbr flux from quasars alone 

The SED of the diffuse MBR resulting from quasars alone 
(black long-dashed line in Fig.^,) is significantly softer 
than a powerlaw. Calculations with such an SED are shown 
in Fig.|21 The calculated AVjiv/^hi ratio assuming 10% 
solar gas lies below the observed value in 0200+015, by a 
factor ~ 10, as shown by the continuous thin line model 
in Fig. |3 Metallicities about solar would be required so 
that the model overlaps the observed column ratio 8 . As for 
0943-242, in absence of definite information about U, the 
absorber's metallicity cannot be constrained any further 
than in the previous powerlaw case covered in Sect. 14.21 

4.3.2. mbr flux from stars and quasars 

In the case in which stars and not just quasars contribute 
to the MBR, the ionizing SED becomes softer (continu- 
ous line in Fig.QJi) and the Ncw/Nm ratio observed in 
0200+015 can now be reproduced using metallicities not 
much above 10% solar, as illustrated by the thin dashed 
line in Fig.QJ The short-long dashed thin line in Fig.0] 
represents the case of an even softer SED in which the 
flux beyond 54 eV has been halved (see dotted line SED 
in Fig.^t). This latter SED hence satisfies our initial goal 
defined in Sect. l4.ll with respect to metallicity. But adopt- 

8 Hence in the thin slab case, the increase in Nciv/Nm pro- 
vided by a quasar-based mbr SED is a factor ten with respect 
to a straight AGN powerlaw. 
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Fig. 3. The column ratio Nqw/Nyu derived from pho- 
toionization by the diffuse mbr due to quasars only, 
as a function of U. The continuous black line and the 
long-dashed black line correspond to the thin absorber 
case assuming metallicities of 0.1 solar and solar, respec- 
tively, while the dotted gray line corresponds to the thick 
absorber case assuming metallicities of 0.01 solar. The 
nomenclature and symbols have the same meaning as in 
Fig. [3 In all figures, black line models only apply to the 
thin absorber objects shown at the top, while gray-line 
models only apply to 0943-242 below. 
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Fig. 4. The column ratio Ncrv/Nm derived from pho- 
toionization by the MBR due to stars and quasars, as a 
function of U. The black short-dashed line corresponds 
to photoionization by the original FGS SED in which the 
mbr flux is generated by quasars as well as stars, as- 
suming a gas metallicity of 0.10 Zq. The black short-long 
dashed line has an SED similar to the previous, but its 
flux beyond 54 eV has been halved (the dotted-line SED 
in Fig.[IJi), owing to a larger contribution by stars. The 
gray dot-dashed line and gray dotted line assume this lat- 
ter SED in the thick slab case, but with metallicities of 0.1 
and 0.01 Zq, respectively. The nomenclature and symbols 
have the same meaning as in Fig.[5] 



ing an SED in which star-forming galaxies are contributing 
more than quasars does not imply that such a distribution 
is typical of the average MBR. It only suggests that this 
SED is valid in the neighborhood of 0200+015. Smette et 
al. (2002) found for instance that the softness of the mbr 
energy distribution (parameter 77, see Sect. l3.2)l presents 
important local variations, with some locations where only 
quasars are apparently contributing while in others there 
appears to be a significant contribution from starbursting 
galaxies. 

In the case of the thicker absorber in 0943-242, there 
is little difference in Nqw/Nhi between the Q+^r sed in 
which the flux beyond 54 eV has been halved and the SED 
produced by quasars only. Compare, for instance, the gray 
dotted lines (0.01 Zq) in FigsElandU 

In summary, a diffuse mbr sustained by quasars and 
star-forming galaxies is quite successful in reproducing the 
observed column ratio in 0200+015 without need for a 
metallicity any higher than ~ 0.10 Zq. 



4.4. Photoionization by local stellar UV 

4.4.1. Photoionization by hot stars with T eff =80 000K 

Using a metal- free stellar atmosphere of 80 000 K and a 
gas metallicity of 4% solar, BG03 obtained a reasonable 
first order fit to the strong lines observed in the unusual 
spectrum of the high redshift LAN. In Fig.|Sl we show that 
the column density ratios of the LAN absorbers can be re- 
produced using a high value of U and an absorption gas 
metallicity of 0.1-0.2 Zq. This range is consistent with 
the comprehensive metallicity determination of the nebu- 
lar emission gas by VM04, that is ~ 10% solar. Given the 
similarity of the LAN column ratio with that of 0200+015, 
we infer that a ATciv/^Vhi ratio of order unity in an hzrg 
is compatible with a stellar SED photoionizing a subso- 
lar metallicity absorber. Hence, the possibility that the 
0200+015 absorber might be photoionized by hot stars 
warrants consideration, since metallicities of only 10% so- 
lar would be needed rather than a value 100 times larger 
favoured by the powerlaw SED (Fig.|2|or J03). The problem 
of stellar continuum detection is discussed in Sect . 15. 2. 21 
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Fig. 5. The column ratio Ncw/Nm derived from pho- 
toionization by a metal-free stellar atmosphere with 
T e ff=80 000K as described in Sect. 13. 51 as a function of U. 
The nomenclature and symbols have the same meaning as 
in Fig.H 



In the case of thicker absorbers, as represented by 
0943-242, the gray line models point to metallicities in the 
range 0.01-0.1 Zq, assuming U < 0.05. Higher values of 
U would imply lower absorber metallicities. Interestingly, 
the stellar (Fig.[5J| and the powerlaw (Fig.|2J models with 
0.01 Z cross the N C w/N m ratio of 0943-242 at a very 
similar U value. This indicates that thicker slabs are much 
less sensitive to the SED's shape. 



4.4.2. Varying the stellar atmosphere temperature 

In Fig.© we explore the effect of varying stellar effective 
temperature. The zero-age metal-free atmospheres used 
correspond to T off of 42 000, 57000, 71000, 80 000 and 
88 000 K (from Schaerer 2002). All models are character- 
ized by a column Nm = 10 148 cm -2 and a gas metallicity 
of 10% solar, appropriate to the LAN. Although temper- 
atures lower than 80 000 K can easily fit the LAN column 
ratio, this would imply too weak HeiiA4686 emission for 
the nebula. At the other temperature end, a T e g as high 
as 88 000 K would require a ten times higher gas metallic- 
ity in order to reproduce the observed column ratio. The 
reason is that, as T e g is increased much beyond 70 000 K, 
the increase in the continuum's hardness causes the slab 
to harbour many ionization stages of carbon (e.g. C +4 
and C +5 ), thereby causing a relative reduction of the C iv 
fraction. Increasing the temperature much beyond 10 5 K 
would result in column ratios approaching those of a pow- 
erlaw. 
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Fig. 6. The column ratio Nciv/Nm derived from pho- 
toionization by stellar energy distributions of varying T e g, 
as a function of U. In all models, the metallicity is 0.1 Zq, 
and the slab opacity 10 14 8 cm -2 . A label shows the T e g 
being considered. The dotted line labelled VM04 corre- 
sponds to photoionization by a stellar cluster of metallic- 
ity 20% solar and 3.4 Myr of age, as described in Sect. 13. 51 
The nomenclature and symbols have the same meaning as 
in Fig.H 

VM04 have modelled the LAN emission line spectrum 
using the ionizing spectrum of an evolved stellar cluster 
in which transient Wolf-Rayet stars can account for the 
nebular He n A4686 line observed in emission. A photoion- 
ization model of the absorber using such an SED is repre- 
sented by the dotted line labelled VM04 in Fig.© The be- 
haviour of the column ratio is similar to that of a 80 000 K 
metal- free star (Fig.JSJ. 

In the case of the hzrg shells, stellar SEDs are possible 
candidates for the ionization of the absorbers (but not of 
their eelr), since they can reproduce the observed column 
ratio using subsolar metallicities. 

5. Estimates of the radii and masses of the shells 
and compatibility with other observables 

Armed with the results of the photoionization calcula- 
tions, we now investigate two of the scenarios in more de- 
tail, namely the cases of ionization by the MBR and by hot 
stars. We focus on their implications for other properties 
of the absorbing shells (e.g. mass, radius and thickness) 
and their compatibility with other observables (e.g. the 
underlying stellar continuum in the case of ionization by 
hot stars, and the strength of the Lya emission). On the 
basis of such considerations we demonstrate that ioniza- 
tion by hot stars is favoured over MBR ionization. Readers 



10 



Binette et al.: hzrg haloes 



who do not wish to follow the argument in full may skip 
over sections 15 . 1 1 and IS~21 and proceed directly to the sum- 
mary in section 15.31 

5.1. The case ofMBR ionization 

We first analyse the possibility of having the MBR ionize 
the haloes. The diffuse mbr is ubiquitous and its intensity 
independent of distance to the hzrg. Therefore, changes 
in excitation (i.e. U) are obtained by varying the gas den- 
sity. Higher shell densities, hence lower U, might explain 
the absence of C iv absorption in many HzRG shells. We 
must thus investigate whether the mbr is strong enough 
to result in an acceptable halo density, because the geo- 
metrical thickness of the shell increases as the density is 
reduced. Below we use such constraints to infer the shell's 
minimum distance from the hzrg and its total mass. 

5.1.1. mbr intensity and shell thickness 

An estimate of the MBR mean intensity is provided by the 
proximity effect, whereby absorbers becoming more ion- 
ized in the vicinity of quasars. We will adopt the value 
J y w 10~ 21 erg cm -2 s _1 Hz -1 sr _1 inferred by Cooke, 
Espey & Carswell (1997) and assume the MBR flux SED 
Q+t't of Fardal etal. (1998), albeit with the flux above 
4Ry divided by two, as studied in Sect . 14. 3 .21 We will 
consider two cases: the weak and the strong absorber 
cases, using 0200+015 and 0943-242 as examples, re- 
spectively. Using the definition 9 of U and mappings Ic 
to integrate the Q+~k SED, we obtain in the optically 
thin case that the total hydrogen density is given by 
n$ in = 2.8 x 10- 4 J_2i{f/o h i in } _1 cm- 3 , where U^ n = 
U thin /0.1 and J_2i = J^/(10- 21 ergcm- 2 s- 1 Hz" 1 sr- 1 ). 
In the case of thicker absorbers with Ah ii 10 18 cm~ 2 , 
because of self-shielding, illumination of a spherical shell 
can only occur from the outside, and the mean intensity 
is approximately half of the previous thin case, such that 
n thick = 2.8 x 10- 3 J_2i{f/o h oo5}~ lcm ~ 3 > where ( for con- 
venience) Uffiftfe = t/ thick /0.005. The photoionization cal- 
culations at constant Am (either 10 14 8 or 10 19 cm~ 2 ) in- 
dicate that the total hydrogen column of the slab can be 
approximated as A H hin ~ 2.1 x 10 19 {Uffi*} 1 - 1 cm" 2 and 
^thick _ 6 3 x 10 20 | f/ thick}i.i cm -2 j reS p ect i ve i y . As ar- 
gued in Sect. 11.31 a shell geometry is more appropriate 
than that of a filled sphere. We therefore introduce an as- 
pect ratio A = Ar/r for the shell, where Ar is the shell 
thickness and r its outer radius, taking the HZRG nucleus 
as the center. Since this ratio is not known, we will de- 
fine an upper limit of A 0.2. Since the total column 
density is given by Nn = rn^A, this limit on A trans- 
lates into a lower limit for the shell radius (i.e. a min- 
imum radius) of > 122 {U^ n } 21 {A , 2 J_ 2 i} _1 kpc 
and rj£ c ck > 365 {U^} 21 {A . 2 kpc, respec- 

9 If the SED was a powerlaw of index a (< 0), we would 
have U ~ — 4-7T J v / (nnhca) — 6.3 10 -5 J_2i/wh, where h is the 
Planck constant and c the speed of light. 



tively, with A a2 = A/0.2. To be definite, we will as- 
sume that both absorbers have the same gas metallicity 
of 0.10 Zq. From Fig.H we read off values of U^ 11 ~ 2 
and ~ 1 for 0200+015 and 0943-242, respectively. 

This translates into minimum radii of 523 and 365 kpc, 
respectively. Hence MBR ionization implies that the shells 
are extremely distant from the background eelr, and this 
appears to be difficult to reconcile with the observations 
which show a distinct transition from sources with radio 
extents < 25 kpc to those with large radio sizes (see sec- 
tion EOJ- 

Since r^ pc oc U 21 , smaller radii follow from assuming 
smaller values of U, which would require that we adopt 
metallicities somewhat larger than 0.1 Zq (larger metal- 
licities shift models to the left in Fig.0}. Uncertainties in 
U (or equivalently in Z) therefore affect our estimates of 
the shell's geometrical thickness significantly. In Sect. l6.ll 
and 16.21 we indicate how detection of the shells in Mgn 
or O vi would help to constrain both Z and U. 

5.1.2. Problem of the large shell masses 

Since the shell masses are given by Mh = 47rr 2 mH-/VH = 
10 _13 r^ JVhMq, assuming they are spherical, we can 
use the previous expressions for the minimum radii to 
derive the following minimum masses Mjj > 3.1 x 
10 10 {[/,5 h 1 in } 5 - 3 {Ao.2J-2i} _2 M Q and M H llick > 8.3 x 
10 12 {U^} 5a {A . 2 J-2i}~ 2 M , respectively. Adopting 
the same estimates of U as above, we derive masses of 
M H hin > 1.2 x 10 12 M Q and M^ ick > 8.3 x 10 12 M Q for 
0200+015 and 0943-242, respectively. At face values, these 
values are excessive and would suggest that MBR ioniza- 
tion is unworkable. On the other hand, given the strong 
dependence of Mh above on poorly determined quantities, 
the upper limits mentioned above are order of magnitude 
estimates and as such do not allow us to completely rule 
out MBR ionization. For instance, reducing both ionization 
parameters by two reduces the 0200+015 and 0943-242 
shell mass estimates to 3.1 x 10 10 and 2.1 x lO n M , 
respectively. Proportionally smaller masses would be im- 
plied if the shells covered only a fraction of 4ir sr. On the 
other hand, if the shells were to be geometrically very thin 
(i.e. A Q . 2 < 0.2), the mass of the 0943-242 shell would be- 
come unreasonably large (> 10 13 Mq). 

In conclusion, the ionization of the shells by the dif- 
fuse mbr would imply that the shells have expanded to 
large distances from the parent hzrg. This would favor 
the 'aging shell' scenario B. A significant problem is that 
the shell mass estimates turn out too large. An alternative 
is that the ionizing radiation is stronger as a result of local 
stellar sources, as discussed below. 

5.2. The case of ionization by local stellar sources 

The similarity of the iVciv /-^hi ratio between the stellar- 
excited LAN nebula and the 0200+015 absorber suggests 
that hot stars could be the ionization source of the hzrgs 
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haloes. Even though the emission-line spectra of the back- 
ground eelr is clearly AGN-like and presumably ionized 
by the hidden quasar, an interesting result of the calcula- 
tions in Sect. 14. 41 is that the column ratios in the two hzrg 
absorbers can be reproduced using a stellar SED and sub- 
solar metallicites, as for the LAN. We now analyse some of 
the implications of this hypothesis. 

5.2.1. Test case: hot stars contributing little to the 

EELR 

The geometry that we envisage is that of a large popu- 
lation of hot stars, possibly distributed uniformly or in 
large aggregates as a result of merging (e.g. the hzrg 
4C41.17; vanBreugel etal. 1997). To simplify the treat- 
ment of the geometrical dilution of the ionizing radiation, 
we will assume that the propagation of the photons is 
approximately radial by the time they reach the interven- 
ing shells. To facilitate the comparison with the previous 
mbr ionization case, we define a reference test case with a 
much higher shell density of 0.01 cm -3 . The photon den- 
sity is set by the relation nn = 10" 2 {£/o.i} _1 c m_3 , which 
is equivalent to having an ionizing flux (reaching the shell) 
36 times higher than that provided by the mbr intensity 
with J-2i = 1, as assumed in Sect. l5.1.D Under the con- 
ditions of this test case, our calculations indicate that the 
ionizing photon flux impinging upon the inner boundary 
of the shells is ipn — 3.0 x 10 7 no. oi?7o.i quanta cm -2 s _1 , 
where n .oi = tih/0.01 represents the shell density. Local 
stellar sources (in contrast to the MBR case) are in better 
accord with the 'inner shell' scenario A, in which the shells 
do not extend further out than about 25kpc in radius, 
i.e. the apparent crossover point between sources with ab- 
sorbers and those without [see J03; W04; Sect . 11.1 1 and the 
superwind-bowshock model of Krause (2005)]. The pho- 
ton luminosity is Airr 2 ^, which can be written as Qh = 
0.224 1 55 r| 5 no. oi Uq.i quanta s _1 , where r25 = rk pc /25. 
The Lya luminosity from recombination alone is given by 
the expression L^ ya = 1.06 x 10~ n e nc b Qn ergs -1 , where 
the conversion factor 10 assumes caseB and a temperature 
of 20 000 K. e ne b is the fraction of photons absorbed and re- 
processed by the emission gas. The leaking fraction 1— e ne b 
for very luminous Hil regions lies in the range 0.3-0.5 
(Beckman etal. 2000 and references therein; Zurita et.al. 
2002; Relaho et al. 2002; Giamanco et al. 2005). To be def- 
inite, we adopt 0.5 and define eg c ^ = e nc b/0.5 to obtain 
that 

L Lya = 0.12 x 10 44 eS c 5 b r 2 5 n .oi U .i ergs" 1 (1) 

for our test case. 

The L^Lya luminosity in the test case should be com- 
pared with the significantly larger eelr L\^ ya luminosi- 
ties of 1.2 x 10 44 and 1.9 x 10 44 ergs -1 , observed in 
0200+015 and 0943-242, respectively. In the case of the 



Under the quoted physical conditions, 65% of the recombi- 
nations lead to Lya photon emission (e.g. Binette etal. 1993). 



LAN, with L Lya — 0.40 x 10 44 ergs -1 , its luminosity 11 
is three times higher than our test case. As for our two 
eelrs, they are brighter in Lya by a factor 5 (0200+015) 
and 300 (0943-242), assuming in Eqn. 1 that Uffi 1 ~ 2 
and f7o!oo5 — 1> respectively. The assumed stellar ioniz- 
ing luminosity is therefore not expected to alter the AGN 
character of the hzrg emission spectrum, even though spe- 
cific emission lines would be subject to a contribution from 
the proposed stellar sources. 

Interestingly, the Lya luminosities of the absorp- 
tion shells themselves are expected to be relatively 
small. We derive Lya luminosities of LLya = 5.5 x 
10 e s heii ^o. oi Lo.ierg s , where e s h n is the fraction of 
ionizing photons absorbed by the shell, a quantity that is 
set by the shell opacity. Our calculations with Ahi of 10 14 8 
and 10 19 cm -2 indicate that e s hcii = 5 x 10 -4 and 0.97 for 
the thin and thick absorbers, respectively. Assuming as in 
Sect. l5~T~Tl that U^ n ~ 2 and U$$g ~ 1 for 0200+015 and 
0943-242, respectively, this translates into luminosities of 
5.5 x 10 40 and 2.7 x 10 41 ergs -1 , respectively. These values 
are negligible with respect to the observed hzrg and LAN 
Lya luminosities. 

As for the masses of the shells, the total column den- 
sities as a function of U in the case of stellar SED are as 
follows: N^ in ~ 1.1 x 10 19 {Uffi*} 1 - 1 cm" 2 and N^ ick ~ 
7.2 x 10 20 {U^} 1 - 1 cm" 2 , assuming the 80000K SED. 
We used the expression M H = 6.25 x 10 -11 r 2 5 A H Mq, 
assuming again that the shells are spherical, to derive 
mass estimates of M$ in = 6.7 x 10 8 r| 5 {U^} 11 M 
and M*[ hick = 4.5 x 10 10 r 2 5 {Utfffife} 1 - 1 M for 0200+015 
and 0943-242, respectively. Hence, the shell masses for the 
test case are quite small in comparison with the mbr case 
because of the smaller radii implied by the stronger ioniz- 
ing flux. By the same token, larger densities are implied, 
which result in shells that are also geometrically very thin 
(A). 2 <C 1). The absorber's density can be quite differ- 
ent than the assumed test case with Uo.oi = 1- The re- 
quired stellar luminosity, however, must then scale in the 
same proportion. For instance, an absorber with density 
0.1cm -3 would require a 10 times higher stellar luminos- 
ity. This would cause the nebular lines to be comparable 



11 The observed I/Lya values quoted above were derived us- 
ing the Lya fluxes reported by V097 and that of Fosbury et 
al. (2003) for the LAN. We assumed the objects to be isotropic 
emitters and corrected the fluxes for Lya absorption due to 
the absorbing shells. For the LAN, L-L ya was divided by 10 
to compensate for the amplification by the gravitational lens. 
We adopted the concordance ACDM cosmology with param- 
eters with fl A = 0.7, Q M = 0.3, h = 0.70 with h = H o /100. 
The isotropic photon luminosities that we infer for 0200+015, 
0943-242 and the LAN are Qn = 1.1 X 10 55 , 1.8 x 10 55 and 
0.37x 10 55 quantas -1 , respectively. These values are lower lim- 
its, since they only represent the fraction absorbed by the gas 
and reprocessed into line emission. To recover the intrinsic Qh, 
they would have to be increased by e nc b -1 , a poorly deter- 
mined quantity in AGN (e n ^ ^ 0.1 : Oke & Korycansky 
1982; Antonucci etal. 1989). 
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in luminosity to the observed eelr which would clearly 
not be desirable. 

5.2.2. On the detection of stellar continuum 

For the case where hot stars alone ionize the foreground 
absorbers, we now estimate the implied stellar flux (or, 
cquivalcntly, the Lya equivalent-width) and compare it 
with the observations, beginning with the LAN. 

Assuming a Salpeter IMF and the SED for an in- 
stantaneous burst of age 3.4 Myr (VM04), we find using 
mappings ic that the rest-frame Lya equivalent-width is 
EWl™ = 190e nc bA. Defining the fraction of ionising 
photons reprocessed by the emission nebula as eg° 5 b = 
Eneb/0.5, we obtain that the (observer-frame) continuum 
flux is F° bs = 0.0105 {eg° 5 b (l + z^y 1 F^A" 1 , where 
F£ya i s the observed line flux, corrected for absorption. 

For the LAN, this implies a 5300 A continuum of 
F° bs = 9.4 x 10 -18 ergcm~ 2 s _1 A -1 , or equivalently, 
8.8 /iJy. (The lens amplification was assumed to be the 
same for both the continuum and the lines). This flux 
is about 30 times higher than the upper limit set by 
Fig. 5 of Fosbury etal. (2003), of w 0.3/xJy. A contin- 
uum was detected at longer wavelengths, but these au- 
thors report that it is consistent with being nebular in 
nature. As a solution, Fosbury et al. (2003) proposed a 
top-heavy IMF. Assuming a single T c g SED of 80 000 K 
(Fig.Hl, we derive EWJ£* = 1335e ncb A, or a continuum 
of 1.3/^Jy(eg c 5 b = 1). Even for a higher T eff of 88 000K, we 
obtain a value of 1.1/iJy, similar to before. The increase 
of a factor 7-8 of the Lya equivalent-widths provided by 
these two SEDs is therefore insufficient. A significantly hot- 
ter stellar SED is therefore required (Fosbury et al. pro- 
posed ~ 10 5 K). Alternative explanations might consist of 
a peculiar dust distribution that preferentially absorbs the 
continuum and thereby increases the observed equivalent- 
width or, as suggested by MV04, differential amplification 
of the lines and the continuum, by the gravitational lens. 
It is interesting to note that the Lya-emitting 'blobs' as- 
sociated with Lyman Break Galaxies likewise do not show 
the expected level for the stellar continuum (Steidel et al. 
2000). For a possible explanation involving significant pop- 
ulations of metal-free stars, see Jimenez & Haiman (2006). 

If we now turn to the two hzrgs, we place upper lim- 
its on the underlying stellar continua by defining lim- 
its on the contribution of hot stars to the eelr Lya 
(which must be much smaller than the AGN contribution). 
Assuming the VM04 SED and the stellar contribution to 
be no more than 10% of the eelr, we derive continuum 
fluxes (at the observed Lya wavelength) of 7.0 x 10~ 19 and 
1.0 x 10- 18 ergcm- 2 s" 1 A" 1 , for 0943-242 and 0200+015, 
respectively 12 . These lie below the upper continuum lim- 

12 It should be emphasized that these are maximum estimates 
of the stellar continuum. We recall that we can let the stel- 
lar continuum be much weaker than the test case explored in 
Sect. l5.2lTl and still have it ionize the shells. Furthermore, if 

U < 0.1, as considered for 0200+015 in Sect . 14.4.11 an even 



its of 3.1 x 10~ 18 and 2.6 x 10~ 18 ergcuT 2 s _1 A -1 , re- 
spectively, as measured by vanOjik (1995). However, us- 
ing recent VLT data from VIMOS-IFU (van Breukelen, 
Jarvis & Venemans 2005), one of us (MJ) reports detec- 
tion of the underlying continuum in 0943-242 at the level 
of 1.7 ± 0.9 x 1Q -I8 ergcm _2 s -1 A -1 , that is, less than 
a factor two above our limit. Vernet etal. (2001) reports 
on the measurement of a far-UV continuum in the form of 
"single peaked sources" . This continuum, however, is 6.6% 
polarized near 1350 A. Vernet etal. (2001) estimates that 
the AGN contributes between 27% and 66% of the con- 
tinuum at 1500 A. After allowing for a 20% contribution 
from the nebular continuum, these authors conclude that 
between 14 and 55% of the unpolarized continuum might 
be due to young stars. The continuum measured by MJ 
is then fortuitously consistent with our upper limit, since 
half of it or less is stellar in origin. We conclude that an 
instantaneous burst with a Salpeter IMF is thus a feasible 
source of ionization for HzRG absorbers. It would in any 
case be difficult to rule it out since continua much weaker 
than assumed in our test case (by a factor ~ 20) would 
still suffice to ionize the absorbers. 

5.3. Summary of current constraints on UzKG halo 
ionization 

To summarise the results of the previous two sections, 
we conclude that ionization by the MBR or by hot stars 
can satisfactorily reproduce the observed Nciv/Nm ra- 
tios without recourse to excessive galaxy-to-galaxy metal- 
licity variations. That was the aim of the photoionization 
modelling as defined in section |4~D 

On closer inspection, however, ionization by the mbr 
leads to excessively large radii for the absorbing shells and, 
by implication, to very large gas masses. This follows be- 
cause the intensity of the MBR, J„, is not a free parameter 
so constraints on U translate directly into constraints on 
halo gas density. Given the latter, the observed require- 
ment for a shell-like geometry translates directly into a 
minimum shell radius from the parent hzrg. For both 
thick and thin absorbers, the minimum radii are of or- 
der several hundred kpc. This is hard to reconcile with 
the observed transition in radio source size between hzrgs 
with and without strong absorption. Scaling as the square 
of the radius, the implied shell masses are also uncomfort- 
ably large. 

The case of ionizing the absorber, but not the eelr, 
by hot stars circumvents the above problems, but at first 
sight raises separate issues of its own. The first is to ensure 
that these hot stars do not overproduce the Lya emission, 
because in hzrgs the eelr is powered by AGN photoion- 
ization or jet interactions. In both 0200+015 and 0943-242 
it was shown that Lya emission from hot stars does not 
significantly contaminate the eelr emission. Potentially 
more serious is the apparent faintness of the stellar con- 
weaker continuum is needed, as shown by the estimates of Lya 
luminosity reported in Sect . 15.2.21 
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ionization parameter U 

Fig. 7. The column ratio A^cii/-/Vhi derived from pho- 
toionization by four different SEDs discussed in Sect.QJ as 
a function of U. The nomenclature and symbols have the 
same meaning as in Fig. [21 

tinuum, which is hard to explain away with peculiar dust 
geometries if the stars ionize gas in the direction of the 
observer. For the LAN Fosbury et al. (2003) appealed to 
hot stars and a top-heavy IMF; for the two HzRGs, con- 
straints on the continuum level below Lya appear to be 
consistent with the levels expected from hot stars. For all 
these reasons, we thus favour hot stars as the more likely 
source of ionization for the hzrg haloes and in the next 
section outline some new diagnostics to test this further. 

6. New diagnostics for future observations 

In order to resolve pending issues such as the size, masses 
and nature of the HzRG haloes, more extensive observa- 
tions are needed and measurements of C IV in absorp- 
tion in other ffiRGs should be attempted. The detection 
of other absorption species would also help to break the 
Z-U degeneracy, as outlined below. 

6.1. Absorption by lower ionization species: Cii\X1335 
and Mgn XX2798 

It would be helpful to detect the absorption of other res- 
onance lines in the spectra of hzrgs, particularly in the 
case of species of lower ionization than C iv. This could be 
used to confirm whether those H I absorbers without C iv 
absorption might simply correspond to shells of lower ion- 
ization (smaller U). Two candidate species are C IIAA1335 
and Mgn AA2798. We report calculations for these two 
resonance lines in Figs [7| and [SJ assuming those SEDs 
that were most successful in reproducing the observed 
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ionization parameter U 

Fig. 8. The column ratio Nugn/Nm derived from pho- 
toionization by four different SEDs discussed in Sect.fJ] as 
a function of U. 

Nciv/Nm ratios. Because both Cll and Mgll are much 
weaker emission lines than Lya, we consider it feasible 
to detect the corresponding absorption doublets only in 
the case of the thicker Hi absorbers. For instance, for an 
absorber with TVhi — 10 19 cm -2 , we expect the Cll and 
Mgn columns to be of order 10 14 and 10 13 cm -2 , respec- 
tively, assuming a metallicity of 0.01 Zq. Interestingly, the 
behaviour of the Ncu/Nhi and Nyigii/Nm ratios is rela- 
tively flat in the strong absorber case, with a dependence 
on U that is much weaker than was the case for C iv. 
This property would facilitate the determination of the 
gas metallicity. A possible strategy would be to use Mg n 
to ascertain the metallicity, and then use the appropriate 
C iv curve to constrain U. 

6.2. Absorption by higher ionization species: 
OviX1035 and NvX1240 

McCarthy (1993) produced a composite optical-UV spec- 
trum of 3CR and 1 Jy sources (redshifts up to 3) that 
is useful for estimating typical strengths of various emis- 
sion lines. Their composite shows that the strongest res- 
onance emission lines in radio-galaxies after Lya and 
CivAA1549 are (in order of decreasing flux) OviA1035, 
O iv+Si IVA1402, N vA1240, Mg n AA2798 and C nAA1335. 
Because O lv+SiivA1402 consists of a blend of two emis- 
sion doublets, it is unlikely that the corresponding ab- 
sorption lines could be disentangled. The other resonance 
lines left to consider are Ovi and Nv. In Fig. 1^1 and 1101 
we present the column ratios Novi/Nm and -/VnvV-^hi, 
respectively, as a function of U. One can see from these fig- 
ures that the ionization parameter could be considerably 
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Fig. 9. The column ratio Novi/Nm derived from pho- 
toionization by four different SEDs discussed in Sect .01 as 
a function of U. 
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ionization parameter U 

Fig. 10. The column ratio N^y/Nm derived from pho- 
toionization by four different SEDs discussed in Sect .01 as 
a function of U. 

better constrained if data on these resonance lines were 
obtained. Thus, obtaining high-resolution optical spectra 
over all emission lines is essential to better constrain the 
properties of these haloes. 
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